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Una prueba empírica de generadores de números 
pseudoaleatorios mediante un proceso de decaimiento 
exponencial. 
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Resumen Las pruebas empíricas que usan procesos o modelos físicos para probar generadores de números 
pseudoaleatorios, complementan las pruebas de aleatoriedad teóricas y han sido usadas con mucho éxito. 
En este trabajo, se presenta una metodología estadística para evaluar la calidad de generadores de números 
pseudoaleatorios, ilustrando el método en el contexto del proceso de decaimiento radiactivo y utilizando 
para ello algunos generadores de uso común en la Física. 

Key words. Aleatorio, pseudoaleatorio, simulación Monte Cario, generador de números pseudoaleatorios, 
Teorema Central del Límite. 
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Empirical tests for pseudorandom number generators based on the use of processes or physical 
models have been successfuly used and are considered as complementary to theoretical test of randomness. 
In this work a statistical methodology for evaluating the quality of pseudorandom number generators is 
presented. The method is illustrated in the context of the so-called exponential decay process, using some 
pseudorandom number generators commonly used in physics. 



PACS. 02.50.-r Probability theory, stochastic processes and statistics ■ 
Monte Cario studies - 05.10.Ln Monte Cario methods. 



02.50.Ng Distribution Theory and 



1. Introducción 

La simulación mediante el Método de Monte Cario [l,2j 
hace uso intensivo de sucesiones de números aleatorios y 
es una técnica estándar ampliamente aplicada desde hace 
mas de 50 años en diversas ramas de la Física, especial- 
mente en la Física de altas energías y en la Física Estadís- 
tica. Además de su aplicabilidad en el método de Monte 
Cario, el uso y obtención de sucesiones de números alea- 
torios constituye hoy en día una investigación muy activa 
con aplicaciones en campos tan diversos como la cripto- 
grafía, integración montecarlo, ecología, identificación bio- 
métrica y aún Inteligencia Artificial [3,4,5,6,7,8,9]. 
Estrictamente hablando, obtener sucesiones de números 
realmente aleatorios implica la utilización de algún fenó- 
meno físico de naturaleza estocástica como el arrojar una 
moneda al aire, el ruido de un circuito electrónico, el de- 
caimiento de un material radioactivo, el conteo de fotones 
mediante detectores centelladores y mas recientemente, se 
han propuesto métodos menos tradicionales basados en 
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fenómenos tales como el flujo turbulento de aire formado 
por el movimiento de los discos duros en una computado- 
ra, péndulos caóticos e incluso del tipo biométricos [TÜ| 
[TT][T2] pero debido a las inherentes dificultades que ofrece 
este enfoque, entre las que podemos mencionar los errores 
sistemáticos introducidos por el arreglo experimental, la 
nula reproducibilidacQ] de la sucesión obtenida, así como 
la baja frecuencia en la generación de números aleatorios, 
han hecho necesaria la búsqueda de otras formas más efi- 
cientes para obtener estos números. 

Desde hace ya algunos años, se utilizan computadoras di- 
gitales para implementar programas a los que llamamos 
generadores de números pseudoaleatorios o simplemente 
generadores, los cuales mediante reglas deterministas y 
operaciones aritméticas muchas veces sencillas, producen 
sucesiones de números que se asemejan en un sentido li- 
mitado ¡13[[14¡. a las obtenidas mediante un experimento 
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1 En la práctica, frecuentemente se requiere usar muchas ve- 
ces una misma sucesión de números aleatorios, lo cual implica 
el tener que almacenar la sucesión completa con todos los in- 
convenientes que esto acarrea. Es en este sentido que usamos 
el término reproducibilidad. 
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aleatorio y que se denominan sucesiones de números pseu- 
doaleatorios. 

Se conocen muchas implementaciones diferentes para ge- 
nerar números pseudoaleatorios [T31Í14] que hacen uso de 
una gran variedad de técnicas y algoritmos que compren- 
den desde el uso de algoritmos de congruencias lineales 
hasta otros asociados con autómatas celulares, algoritmos 
de criptografía de curvas elípticas, etcétera [T5lll6[ fT7]. 
La actual utilización de series muy grandes de números 
pseudoaleatorios en muchas aplicaciones, así como algu- 
nos episodios de resultados dudosos, obtenidos debido a 
la baja calidad de los generadores utilizados [T4lfT8lfT9l 
120] . ha fortalecido la necesidad de contar con mejores y 
cada vez más eficientes pruebas de la calidad. El campo 
de investigación de las pruebas de calidad de generadores 
de números pseudoaleatorios (y por supuesto, también de 
su implementación) , es tan activo que prácticamente no 
hay mes en el que no se reporten en la literatura científica 
nuevas pruebas de calidad que utilizan una gran variedad 
de criterios y técnicas (teoría de la información, técnicas 
estadísticas, power spectrum, gambling tests, sistemas fí- 
sicos, entropía, etc) [3ll2ip22p23p 24 p 25 p26p27p8] . 

Las pruebas de calidad de los generadores de números 
pseudoaleatorios se pueden dividir en: 

■ Pruebas teóricas: Se realizan estudiando los algoritmos 
generadores de números pseudoaleatorios mediante el 
uso de herramientas como la teoría de números. Es- 
tos tipos de pruebas son útiles por su generalidad y 
están basadas en el estudio de algunas propiedades ta- 
les como la longitud del periodo de la secuencia y la 
uniformidad del algoritmo. 

■ Pruebas empíricas: Estas pruebas se concentran en las 
sucesiones de números pseudoaleatorios y sus propie- 
dades. Son usadas para encontrar correlaciones loca- 
les no triviales presentes en las sucesiones de números 
pseudoaleatorios y mostrar aspectos desapercibidos en 
las pruebas teóricas. 

En los últimos años se han explorado pruebas empíricas 
basadas en procesos físicos [|29,30J, principalmente la ca- 
minata aleatoria y el modelo de Ising [25,26,29,3Ó"ll31|. 

En este trabajo, se presenta un criterio para evaluar a un 
generador de números pseudoaleatorios dado, sugiriendo 
que su calidad está directamente relacionada con la ca- 
pacidad del mismo para producir simulaciones que repro- 
duzcan el comportamiento y propiedades teóricas de un 
proceso o modelo de referencia. En otras palabras, nues- 
tro enfoque radica en valorar la calidad de un generador 
mediante una medida única y global de congruencia entre 
las propiedades teóricas del modelo y las producidas por 
el generador, en lugar de aplicar un conjunto de pruebas 
de aleatoriedad aisladas que se basan en desviaciones par- 
ticulares. Esta idea es aplicada, a manera de ilustración, 
para evaluar la calidad de algunos generadores con respec- 
to a su capacidad de simular un proceso de decaimiento 
radiactivo. En la siguiente sección haremos un breve repa- 
so de la ley del decaimiento radioactivo y de su simulación 



mediante el Método de Monte Cario; en la sección [3] pre- 
sentamos los generadores que se probarán en el presente 
trabajo; en la Sección |4] explicamos los criterios de nuestro 
método de prueba de generadores de números pseudoalea- 
torios y en la sección [5] ilustramos la aplicación de nuestro 
método. Finalmente, en la sección [6] se hace una discusión 
de los resultados obtenidos. 

2. La Ley del Decaimiento Radioactivo 

Por razones de autocontención y con el fin de establecer 
la notación usada en este trabajo, revisemos muy breve- 
mente [32j la ley de decaimiento radioactivo: 
Considérese al tiempo t = 0, una muestra grande con A^O) 
partículas inestables y supongamos que cada partícula tie- 
ne una probabilidad XAt de decaer durante el intervalo de 
tiempo pequeño At. Aquí A es la constante llamada razón 
de decaimiento. Si tenemos un número N(t) de partículas 
sin decaer al tiempo t, entonces XNAt partículas decae- 
rán durante el intervalo de tiempo [t, t +At\. Claramente, 
esto decrecería el número de partículas que aún no decaen 
en: 

AN = -XNAt (1) 

Cuando el intervalo de observación en el tiempo At tiende 
a cero, podemos integrar la ecuación diferencial resultante 
obteniendo la ley de decaimiento exponencial: 



N(t) = N(0)e 
Las unidades de A son tiempo -1 . 



-At 



(2) 



Se define la semi-vida como el intervalo de tiempo en el 
cual la mitad de la muestra inicial habrá decaído: 

N(T 1/2 ) = N(0)/2 , es decir: 



T 1/2 = — = 



ln2 0,693 



A 



(3) 



Finalmente, otra definición importante a considerar es la 
de la vida media r. Esta se define como el recíproco de la 
razón de decaimiento: 



1 

A 



(4) 



2.1. Algoritmo para Simular la Ley de Decamiento 
Radiactivo 

El algoritmo utilizado para llevar a cabo la simulación del 
decaimiento radioactivo es el siguiente: 

1. Se eligen los valores iniciales del decaimiento (cantidad 
de partículas iniciales, valor de A, el intervalo de tiem- 
po que se simulará [0,T] ) y el tamaño del incremento 
en el tiempo At. 

2. Para cada valor discreto del tiempo i¿e[0, T], i = 1,2, ..n; 
considérense todas las N*(ti) partículas que no han de- 
caído hasta ese instante. Por medio del generador de 
números aleatorios bajo prueba, genérese un número 
pseudoaleatorio. 
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3. Si el número pseudoaleatorio generado en el paso an- 
terior es menor que lambda, la partícula decae, y que- 
darán iV*(í¿-|_i) = N*(ti) — 1 partículas al tiempo í¿+i. 

4. Repítase el procedimiento a partir del paso 2 para to- 
das las demás N*(U) — 1 partículas y para cada í¿ en 
el intervalo especificado. Termínese el procedimiento 
cuando el número de partículas restantes sea menor 
que 10. (Esto con el fin de evitar fluctuaciones estadís- 
ticas demasiado grandes). 

El símbolo * es utilizado de aquí en adelante para indicar 
que los anteriores valores son los valores producidos por la 
simulación y no por la ley (J2¡) - En la figura Q] se muestra 
el resultado de realizar esta simulación 10 000 veces, con 
razón de decaimiento A = 0.035 y con un número inicial 
de partículas de Nq = 1000. 

N* 



1000 E" 




Tiempo 

Figura 1. 10 000 simulaciones del decaimiento radioactivo con 
A = 0.035 , iVo = 1000 partículas y At = 0.001; la figura con- 
tiene 2 x 10 9 puntos y se obtuvo usando el generador Random3 
descrito en la siguiente sección. El panel superior muestra la 
distribución de N*-N(t) para 1,25, 50 y 75 u. de tiempo. Para 
un buen generador y como consecuencia del Teorema Central 
del Límite, estos valores se distribuirán gausianamente, lo cual 
es uno de los criterios de nuestra prueba, como se explicará 
más adelante en la sección [4] 



3. Generadores 

En este trabajo, y para ilustrar nuestro método se utili- 
zaron los generadores ranO, ran48 y los generadores im- 
plementados en las clases basadas en C++ TRandom, 
TRandom2 y TRandom3 y que son específicos del am- 
biente ROOT [39J, de amplia utilización en física de altas 
energías. Por simplicidad los llamaremos Randoml, Ran- 
dom2 y Random3 y los explicaremos a continuación: 

■ ran48: Pertenece a la familia de funciones rand48() 
generadoras de números aleatorios del lenguaje C es- 
tándar y es un generador MLCG (Multiplicative Li- 
near Congruential) , es decir, usa la relación recursi- 
va r¿ = (a + c)mod m, para ¿=1,2,3... Donde ro, 
a, m y c son números enteros de 48 bits con valores 
a = 25214903917, c = 11 y m = 2 48 . 



■ ranO: Es un generador cuya implementación actualiza- 
da puede encontrarse en [14] y que fue propuesto por 
Park & Miller [33|I34| como un modelo de generador 
con calidad mínima. La implementación usada de este 
algoritmo genera números pseudoaleatorios con distri- 
bución uniforme entre y 1. Se trata de un algoritmo 
lineal congruencial, donde, siguiendo la notación de la 
entrada anterior se hace a = 7 5 = 16807, c = y 
m = 2 31 - 1 = 2147483647. 

■ Randoml: La clase TRandoml, implementa la función 
Rndm(), que es una traducción y actualización a C++ 
de la función fortran generadora de pseudoaleatorios 
RNDM, la cual pertenece a las famosas bibliotecas 
científicas CERNLIB [35j. Randoml genera números 
pseudoaleatorios mediante el método de congruencias 
lineales, con una distribución uniforme en el interva- 
lo [0,1] con un periodo de 10 8 . Los parámetros usados 
por esta función son a — 5 15 = 30517578125, c = y 
m = 2 47 = 140737488355328. 

■ Random2: Implementa un generador de números pseu- 
doaleatorios basado en el el método de Tausworthe 
máximamente equidistribuido por L'Ecuyer y usa 3 
palabras por estado. Toda la información sobre este 
generador se puede consultar en [36p7]. El periodo de 
Random2 es 2 88 = 309485009821345068724781056. 

■ El generador Random3, se basa en el método "Mersen- 
ne Twister" en la implementación de M. Matsumoto 
y T. Nishimura, [38J. Random3 es un generador 623- 
dimensional que produce números equiprobables en el 
intervalo [0,1]. Su gran ventaja es su periodo largo 

, aunque su desventaja es su estado interno 
relativamente grande de 624 enteros. 

Los códigos fuente de los tres últimos generadores se pue- 
den consultar en la página de ROOT[39j. 

3.1. Breve estudio preliminar de la calidad de los 
generadores usados en este trabajo 

Con el objetivo de hacer una comparación independiente 
de nuestra prueba, el cuadro Q] muestra los resultados de 
aplicar 5 sencillas pruebas de aleatoriedad a los genera- 
dores usados en nuestro estudio. Estos se realizaron sobre 
una pequeña muestra de 30 000 eventos para cada gene- 
rador y consistieron en: 

■ Cálculo de ¡x, el valor esperado promedio. Óptimo va- 
lor /i = 0.5. 

■ Cálculo de la desviación estándar. Óptimo valor: 
a = yp^oi 0.28861. 

■ Prueba de Anderson-Darling [40]. El mejor generador 
es aquel que minimize al estadístico A 2 , el cual nos da 
la distancia entre la función de distribución empírica 
y la teórica, en este caso la uniforme en (0,1). 
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■ Estimación de n = 3.14159.... 

■ Cálculo del coeficiente de correlación de Pearson C p . 
Óptimo valor C p = 0. 

■ Compresión de los archivos de la muestras usando el 
algoritmo de Lempel & Ziv [41] y comparando el gra- 
do de compresión con una muestra del mismo tamaño 
de números verdaderamente aleatorios obtenidas de la 
medición del ruido atmosférico ¡42] . Óptimo (ideal) va- 
lor 0%. 



Generador 


/i ± 0,002 


o- ±0,0012 


Á 2 


Randoml 


0.500 


0.2885 


1.145 


Random2 


0.502 


0.2889 


0.280 


Random3 


0.496 


0.2891 


0.523 


drand48 


0.502 


0.2877 


0.859 


RanO 


0.500 


0.2885 


1.048 



Generador 


7T 


Cp 


Compresión 


Randoml 


3.1443 


0.0019 


1% 


Random2 


3.1404 


0.0003 


1% 


Random3 


3.1364 


0.0011 


1% 


drand48 


3.1308 


0.0018 


1% 


RanO 


3.1343 


0.0089 


1% 



Cuadro 1. Algunas sencillas pruebas de calidad realizadas a 
una muestra de cada generador de 30 000 números aleatorios 
uniformemente distribuidos en el intervalo (0,1). Cuadro su- 
perior: Test de uniformidad de los generadores. 
Cuadro inferior: Cálculo de ty mediante el método de Monte 
Cario, Coeficiente de autocorrelación y compresibilidad de la 
muestras generadas. Los parámetros se explican en el final de 
la presente sección. 



4. Criterios de calidad 

Como se indicó en la sección 1, el criterio básico aquí uti- 
lizado para valorar la calidad de un generador en el con- 
texto de un modelo dado, se relaciona directamente con 
la capacidad del generador para reproducirlo y es en ese 
sentido que presentamos una medida estadística de con- 
gruencia, fundamentada en las propiedades teóricas del 
proceso bajo simulación. En nuestro caso, consideraremos 
tres características: 

1. Normalidad conjunta de los valores simulados. 

2. Convergencia de los promedios de los valores simula- 
dos a sus promedios teóricos (convergencia de primeros 
momentos) . 

3. Convergencia de las varianzas y covarianzas muéstrales 
a sus valores teóricos (convergencia de segundos mo- 
mentos). 



Bajo un proceso de decaimiento exponencial como el que 
aquí se ha referido, puede verificarse que el valor esperado 
y la varianza de N* (ti), están dados por 

H(U) = E [N* (ti)} = N(0) exp(-Ati), 

U 

o 2 (ti) = V [N* (U)} = N(0)X ¿ exp [-A (U + j - 1)] , 

3=1 

Además, la covarianza entre las dos variables aleatorias 
N*(ti) y N*(tj) se expresa mediante: 

(j(U,t 3 ) = (l-\) t¡ - u o 2 (U), 

donde í¿ < tj y hij = tj — U, i, j = 1, 2, n, con n igual 
a nuestro número de observaciones. 

Para cada uno de los n instantes de tiempo t%, . . . ,t n en el 
intervalo [0, T], denotemos por N* = [N* (í x ) , . . . ,N* (t n )}' 
al vector aleatorio que representa los valores del proceso. 
El vector N* tiene entonces como valor esperado, al vec- 
tor ¡jl = [¡i (ti) , . . . , // (í„)] y por matriz de covarianzas 
T=[a(t h t j )}. 

Bajo el supuesto de normalidad conjunta, la cantidad 

D = [N*- M ]'r- 1 [N*- M ], 
tiene una distribución y 2 con n grados de libertad. 

Dadas ra realizaciones del proceso, denotemos por Ñ* al 
vector de promedios 

Ñ* = [Ñ(h),...,Ñ(t n )]' , 

m 

donde N* (U) = i ¿ Nj (U) y N* (U) denota el j-ésimo 

3=1 

valor simulado en el punto í¿. 

Independientemente del supuesto de normalidad del vec- 
tor N*, por aplicación del Teorema Central del Límite, el 
vector y/m (Ñ* — /i) converge en distribución a la de una 
variable aleatoria normal n— variada con vector de medias 
y matriz de covarianzas T; de ahí que la cantidad 

D OT = m(Ñ*-/x) / r- 1 (Ñ*- M ), 

tiene como distribución límite, una x 2 con n grados de 
libertad. 

D m es una medida estadística del ajuste entre los valo- 
res teóricos de los parámetros de primero y segundo orden 
derivados del modelo, en este caso el de decaimiento expo- 
nencial, y sus contrapartes muéstrales obtenidos a partir 
de la simulación. Esta medida es la que en este trabajo se 
propone como criterio cuantitativo para la comparación de 
los generadores, ya que su valor se incrementa en presencia 
de desviaciones en los valores teóricos de los parámetros o 
bien en presencia de desviaciones de su estructura teórica 
de covarianzas. En el contexto de nuestro modelo, el mejor 
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generador será aquél que produzca el menor valor de D m . 
Cabe señalar que, con este enfoque, la(s) causa(s) específi- 
ca^) de las desviaciones que pudieran presentarse debidas 
a deficiencias del generador (como correlaciones, ciclos o 
cortos períodos, etc) nos resultan irrelevantes en tanto que 
nuestro interés radica en detectar su efecto final, observa- 
do como incongruencia estadística de las simulaciones con 
respecto a las propiedades del proceso que se simula. 

5. Implementación del Método 

Ilustremos ahora el procedimiento anteriormente descrito, 
considerando un proceso con los siguientes parámetros: 
N(0) = 1000 partículas iniciales, constante de decaimiento 
A=0.035, Al = 0.0025 y m = 2000 simulaciones para cada 
generador. Los n — 6 instantes observados en el intervalo 
[0, 250] fueron h = 1, t 2 = 25, t 3 = 50, U = 75, t 5 = 100 
y í 6 = 125. 

El cuadro [2] muestra los valores calculados de D m para 
cada uno de los generadores y sus valores de probabilidad 
(p) asociados, con base en la distribución % 2 con 6 grados 
de libertad. 



Generador 


D m 


Probabilidad 


drand48 


2.914 


0.820 


ranO 


132.79 


0.00 


Randoml 


14.500 


0.025 


Random2 


2.634 


0.85 


Random3 


1.223 


0.976 



Cuadro 2. Valores calculados de la estadística D„ 
uno de los generadores examinados. 



para cada 



Para estos valores, el vector de medias y la matriz de cova- 
rianzas teóricos del modelo, son: 

fi = [965,605 416,862 173,774 72,440 30,197 12,588] 



33,796 14,372 5,898 2,420 0,9933 0,4076 
247,367 101,514 41,659 17,096 7,016 
146,104 59,958 24,605 10,097 
68,375 28,060 11,515 
29,801 12,230 
12,649 



r 



El menor valor de D m corresponde al calculado para el 
generador Random3, por lo que concluimos que éste re- 
produce el modelo de decaimiento radiactivo con mayor 
precisión. Aunque el generador Random2 presentó un va- 
lor mayor de D m , el valor de probabilidad asociado nos 
indica que, aproximadamente un 85 % de las veces en que 
se simule este proceso bajo las mismas condiciones, po- 
dríamos esperar un valor mayor o igual a 2.634, por lo que 
el valor obtenido es suficientemente pequeño para concluir 
que la simulación es muy aceptable. Por el contrario, RanO 
y Randoml son los peores generadores según nuestro cri- 
terio, el resultado producido por el generador Randoml, 
nos indica que la probabilidad de obtener un valor tan 



grande como D m =14.500 es del orden de 1 en 40, signi- 
ficando con esto que un valor tan improbable nos lleva 
necesariamente, a concluir que este generador no simula 
satisfactoriamente nuestro proceso. RanO es aún peor con 
un valor de D m casi un orden de magnitud mas grande 
que el de Randoml y con probabilidad cero. 
Las figuras[2]a[6]muestran los histogramas de los valores D 
bajo las mismas condiciones de simulación para cada uno 
de los generadores bajo estudio, recordando que para cada 
generador, éstos corresponden a 2000 simulaciones del de- 
caimiento radioactivo con A = 0.035, un número de 1000 
partículas iniciales e intervalo de tiempo de Al= 0.0025 
unidades de tiempo. En ellos se aprecia una congruencia 
del valor obtenido de D m con el ajuste X( 6 ) ilustrado por 
la linea continua en cada una de ellas. 
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Dm = 2.914 
P = 0.820 
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Figura 2. Histograma de los valores de D para drand48. Se 
muestra el resultado de 2000 simulaciones. La linea sólida re- 
presenta la densidad X^6)> Q ue teóricamente tiene D. En esta 
prueba obtenemos los valores D m = 2.914 con probabilidad p 
= 0.820. 
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Figura 3. Histograma de D para ranO. La linea sólida repre- 
senta la densidad X(6); Q ue teóricamente tiene D; note que el 
ajuste no parece satisfactorio, lo cual es claro de los valores ob- 
tenidos de Dm = 132.79 y probabilidad p = 0.0 en esta prueba. 

Es importante destacar que, por el Teorema Central del 
Límite, la distribución de D m es poco sensible a desviacio- 
nes de la normalidad de los vectores N*(í) y consecuente- 
mente a la distribución original de D; sin embargo, a partir 
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Figura 4. Histograma de D para el generador Randoml. 
Como antes, la linea sólida representa la densidad xfa) teó- 
rica de D- note que el ajuste no es satisfactorio, lo cual puede 
verse de valores obtenidos en esta prueba, los cuales son D m 
= 14.500 y una probabilidad p = 0.025. 
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Figura 5. Histograma de D para el el generador Random2 y 
su comparación con la linea sólida que representa la densidad 
de X(6) j ci ue teóricamente tiene D; en este caso, el ajuste parece 
adecuado y obtuvimos los valores de D m = 2.634 y p = 0.85 
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Figura 6. Histograma de D corresponediente al generador 
Random3, con A = 0.035 y n =6. , que teóricamente tiene D; 
para este generador, el ajuste de x%) es e ' mejor de todos, lo 
cual se ve númericamente de los valores D m = 1.223 y p = 
0.976. 



de nuestros resultados, podría percibirse incorrectamente 
que D m es sensible a desviaciones en la distribución teóri- 
ca de los valores D, por lo que es necesario aclarar que la 
congruencia observada entre nuestra medida de ajuste D m 
y las variaciones en la distribución de D con respecto a la 
distribución xfe) > como la que se aprecia en la figura[4j son 
consecuencia de diferencias entre los valores teóricos de los 
parámetros y los producidos por la simulaciones, mismas 
que se ven reflejadas como irregularidades en el histogra- 
ma respectivo. Es, en este sentido específico, que es posible 
concluir que D rn mide indirectamente esas irregularidades 
en la distribución de D, condensando el comportamiento 
general en una sola cantidad. 



6. Conclusiones 

En este trabajo se propone y se fundamenta teóricamente 
una prueba empírica de generadores de números pseudo- 
aleatorios basada en el proceso estocástico de decaimiento 
exponencial. Específicamente hacemos uso de la distribu- 
ción del número de partículas sobrevivientes y su norma- 
lidad en corcondancia con el teorema central del límite. 
Nuestra prueba, tal como debiera esperarse para este tipo 
de criterios de aleatoriedad, posee la característica de ser 
sensible tanto a desviaciones distribucionales como a las 
de los parámetros del modelo usados en la simulación. 
Se ilustra este método mediante su aplicación a cinco gene- 
radores conocidos: tres generadores congruenciales lineales 
con distintos parámetros, un Tausworthe y un Mersenne 
Twister, siendo el mejor de todos según los criterios de 
nuestra prueba este último. 

Nuestro método no resulta de difícil implementación, solo 
hay que tener cuidado al realizar las simulaciones comple- 
tas del proceso de decaimiento radioactivo, con las dificul- 
tades computacionales que esto implica, especialmente con 
lo que respecta al tamaño del intervalo de tiempo elegido 
(véase [43] pg. 53). Por otro lado, es posible, consideran- 
do la longitud del periodo de los generadores modernos 
usar una muestra relativamente pequeña de sus valores 
En este trabajo, para obtener la muestra a analizar con 
nuestro método, se hicieron 2000 simulaciones del decai- 
miento radioactivo para cada generador, esto es en total 
se simularon 50 000 000 eventos y de estos seleccionamos 
2000 valores para cada uno de los seis tiempos distintos 
(1,25,50,75,100 y 125 unidades recíprocas de A). Es de- 
cir, en nuestra prueba y para cada generador, de los 50 
000 000 eventos simulados un total 12 000 eventos fueron 
analizados. 
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